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Abstract 

We reexamine the top quark production cross section at the Tevatron and 
LHC, in the hght of recent progress on the resummation of logarithmic soft 
gluon corrections. We find that resummation effects are much smaller than 
previously thought. We also compute Coulombic threshold effects, and find 
them negligible. We update the discussion of uncertainties due to scale de- 
pendence, the value of the strong coupling constant, and the parton density 
parametrization. Our current best estimate of the top production cross section 
at the Tevatron and its error is a{mt = 175 GeV) = 4.75]'^g'g2- 
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1. Introduction 



We present in this letter a theoretical reassessment of the evaluation of the top 
quark production cross section in high-energy hadronic collisions. The top quark 
having been found [jl|, ||, the comparison between its observed production properties 
and those expected from the Standard Model will be an important probe of the pos- 
sible existence of new phenomena. One of the most important tests to be performed 
concerns the total production cross section. This is the most inclusive quantity avail- 
able, and is a priori the least sensitive to a detailed understanding of the higher 
order corrections which influence the evolution of the initial and final states. Within 
QCD, one expects the perturbative expansion in powers of as{mt) to be well behaved 
and to provide an accurate estimate of the total cross section already at low orders. 
In particular, the first estimates of the total production cross section using the full 
next-to-leading-order (NLO) matrix elements [§, H, ^ gave an increase (relative to the 
Born result) of the order of 30% for masses above 100 GeV. The residual perturbative 
QCD uncertainty, evaluated by varying the renormalization and factorization scales, 
was shown to be no larger than 10%. The choice of parametrization for the input 
parton densities was also shown to give effects of this order of magnitude, by using 
the available sets. 

It was pointed out in ref. that logarithmic contributions associated to the 
emission of soft gluons from the initial state could significantly enhance the NLO 
result. An independent study of the soft gluon resummation in top production has 
appeared recently . These studies are based on the resummation formulae obtained 
in the pioneering works of refs. ||^, |^, . Using these formulae for explicit calculations 
requires the choice of a prescription in order to bypass the Landau-pole singularity 
which is exposed by the integration of the QCD running coupling over gluon energies 
of the order of Aqcd- In particular it was shown in ref. that, within the proposed 
formalism, one can get sensible results only by cutting off the soft gluon emission at 
scales below 0.1 0.2 rrit. The choice made in ref. follows instead a suggestion 



put forward in ref. [jTT| . According to this prescription the Landau pole appearing in 
the Mellin transform of the resummed hard cross section should be integrated over 
in a principal value sense. It was argued in ref. 0] that this prescription results in 
non-perturbative power-suppressed ambiguities scaling like Aqcd/Q in the hard scale 
Q typical of the process under consideration. The net effect of gluon resummation 
evaluated in ref. amounts to an increase by approximately 10% of the NLO cross 
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section. 



In a companion paper [|I2] we will show that the approaches of refs. overesti- 
mate both the gluon resummation contribution and the associated non-perturbative 
residual uncertainty. This is a consequence of the fact that their resummation for- 
mulae introduce unjustified factorially growing terms in the perturbative expansion. 



In ref. |12| we will also show that a more natural prescription exists, in which these 
terms are not present. In the following, we will briefly summarize the basis of our 
criticism and the definition and main properties of the new proposed resummation 
prescription. We will then present a numerical study indicating that the impact of 
resummation on the total top production cross section is of the order of a percent, 
much smaller than previously thought. We will also show that threshold corrections 
due to higher order Coulomb effects are likewise negligible. We will conclude this 
study with an updated analysis of the current theoretical uncertainties coming from 
the scale dependence, from the choice of partonic densities and from the value of as ■ 

It should be pointed out that this study is performed within the strict domain of 
the Standard Model. Corrections to the top cross section much larger than the Stan- 
dard Model QCD uncertainties can be obtained in the presence of new phenomena. 



For a partial list of specific examples, see e.g. refs. [13 



2. The minimal prescription for the resummation 

The heavy flavour production cross section in the Born approximation is given by 
the formula 

a= E f dx,dx2F^''Hx,)F^''^\x2) (-^) , p = ^ , (2.1) 

where m is the mass of the heavy quark, and S is the square of the total centre-of- 
mass energy. The functions are the parton densities for parton k in hadron 

i, and are evaluated at a scale fi of the order of the heavy quark mass m. Explicit 
formulae for the partonic cross section a can be found for example in ref. 0. The 
corresponding space formula is 

^= E l^j'^^^ F^Ntl Fj^tl a^;/^ p-^ dN , (2.2) 

— ZTTt JC—ioo 
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where the Mellin transforms are defined in the following way 



Jo X 



dz 



z 



z^'aiz) 



The integration contour lies to the right of all singularities. 
We use the following resummed formula 



a 



(res) 



E 



where, in the MS scheme 

, (y) 
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(2.3) 



(2.4) 



(2.5) 

(2.6) 
(2.7) 

(2.8) 



In formula ( p.5|) the strong coupling constant is evaluated at a scale /x of the order of 
the heavy quark mass m. The factor In A^^'^'* resums all the leading logarithmic terms 
In*^^^ due to soft gluon emission. The contour in eq. (|2.4| ) is chosen between 
the cuts on the negative N axis and the singularity at hoaglogN = 1/2, which is a 
Landau pole. We will call this the minimal prescription (MP). In ref. [0 we will 
show that this prescription enjoys the following remarkable properties: 



• It is accurate at the leading log level in the threshold limit. 

• If we expand the MP formula in powers of a(m), the expansion is asymptotic 
to the full formula. 

• The power expansion of the MP formula does not have any factorially growing 
coefficients, and therefore is free of ambiguities of the order of powers of A/m. 
The ambiguity associated with the asymptotic expansion of the MP resumma- 
tion is exponentially suppressed, being of the order e~^^^~P^^. 
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We claim that previous calculations of resummation effects have a corresponding 
power expansion which does have factorially growing terms. These terms are spurious. 
They arise either from certain approximations performed when going from the N space 
to the X space formulae, or from the integration of the running coupling constant down 
to the Landau singularity. 



Our result is consistent with the findings of ref. [14], indicating that no power-like 
1/Q ambiguity emerges from the resummation of multiple gluon emission in Drell-Yan 
processes. 

A full discussion of these issues will appear in a companion paper |12|, where we 
will also study the implications of this prescription for Drell-Yan, heavy quark and 
jet production in hadronic collisions. 

In the present work we will report on results obtained with our MP formula of 
eq. (|2.4|) , which was implemented without any further approximation in a numerical 
program. 



3. Results 

In this section we present the results obtained for a very heavy quark, and in 
particular for the top cross section. First of all, we remind the reader that large soft 
gluon effects are present not only in the production cross section, but also in the 
deep-inelastic processes that are used to determine the structure functions. In order 
to perform a complete study of the resummation one should therefore use resummed 
formulae also when fitting deep-inelastic scattering data, direct photon production 
data, and in general all phenomena which are used to constrain the parton densities. 
In the present work we will assume that the structure functions have been properly 
extracted from data (including the resummation effects) and assess the significance of 
the resummation for heavy flavour production. Effects of comparable size may arise 
from refitting the structure functions using resummed formulae. Our conclusions will 
mostly be based upon the fact that the resummation effects we find are in fact very 
small, and can be neglected. 

The importance of the resummation effects is illustrated in fig. |1| where we plot 
the quantities 

^gg ^qq (^gg + Sqq , , 

Jga) ' (qq) ' (99) ^"^'^^ 
'-'nlo "^NLO "^NLO "t "^nlo 
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Figure 1: Contribution of gluon resummation at order and higher, rela- 
tive to the NLO result, for the individual subprocesses and for the total, as a 
function of the top mass in pp collisions at 1.8 TeV. 



where 5 is equal to our MP resummed hadronic cross section ( |2.4| ) in which the terms 
of order a] and have been subtracted, and cr(NLO) is the full hadronic NLO cross 
section [Q]. In other words, what we display is the relative contribution to the total 
cross section coming from corrections of order higher than those accounted for by the 
full NLO result. We show the results for the qq and gg production channels separately, 
as well as for the total. The results are shown at = L8 TeV function of the 
top mass in the range 100 < mt(GeV) < 600. The wide mass range is chosen for the 
sole purpose of illustration. We chose a common renormalization and factorization 



scale n = rrit, and parton densities from the CTEQl set |]T5 



The plot clearly indicates that the contribution of soft-gluon resummation is neg- 
ligible, unless the top mass approaches the total hadronic energy. This result is 
consistent with the expectation that the soft gluon enhancement of the production 
cross sections should be relevant only very close to threshold. For a top mass of 
175 GeV the increase due to this effect is below 1%. We also note that the effect 
of resummation is potentially much larger for the gg channel, because of the Cx/C-p 



enhancement in the exponent of the resummation formula (see eq. |2.6| ). It turns out, 
however, that when the gg subprocess becomes important [i.e. for small top masses) 
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Figure 2: Contribution of gluon resummation at order aj. and higher, relative 
to the NLO result, for the individual channels and for the total, as a function 
of the CM energy in pp collisions. 



we are so far from threshold that the resummation effect is very small. Conversely, 
for large masses, the quark component, which has smaller resummation effects, dom- 
inates. 

In fig. |]we show the same quantities, at fixed top mass, versus the CM energy. We 
see that at the LHC, where production is dominated by the gg initial state, we are far 
enough from the hadronic threshold, so that the resummation effects are again small. 
We conclude therefore that the effects of gluon resummation can be neglected at the 
present stage, as they are much smaller than the uncertainties present in the NLO 
determination. Note that, due to the dominance of the gg subprocess, the relative 
importance of the resummation effects at =1.8 TeV is larger in pp collisions than 
in pp. 

Beside the soft gluon emission effects, also Coulomb effects may enhance or deplete 
the cross section near threshold |jT6|, We have calculated these effects in the 
following way. We have separated the partonic Born cross section formulae into their 
colour singlet and octet components 

^(9s) = afgf + (3.2) 
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'(8) 

where a'^^^'Sg) be found in refs. Il3|, and 



cr(««") = afsf (3.3) 



(1) ~ m2 384 

where (3 = a/1 — p. The Coulomb-resummed cross section is given as 



(3.4) 



(J (p)=cr(8)(p) 7 77—577 - + a(i)(pj- — . (3.5) 

exp(7ras/(6p)) — 1 1 — exp(— 47ras/(3p)) 



We do not include bound state effects, which, as shown in ref. |T^, are much smaller. 
In fig. ^ we have plotted the hadronic quantity 

^Coul 



Cr(NLO) 



(3.6) 



which (similarly to the case of soft gluon effects) represents the relative correction due 
to the resummation of Coulomb effects not already included into the NLO results. 
We see that these effects are similar in magnitude to the soft-gluon effects, and fully 
negligible at the Tevatron. Since at the LHC we are further away from the threshold 
region, we conclude that also there they will be negligible. Therefore, the following 
estimates of the top cross section will not include neither soft gluon effects, nor 
Coulombic ones. Electroweak corrections to top production in hadronic collisions 



have been considered in refs. [|l8l. They range from -0.97% to -1.74% of the born 



cross section, for a higgs mass of 60 and 1000 GeV respectively [jT9|. They will also 
not be included in the following estimates. 

We now proceed to review the current theoretical uncertainties that arise at NLO. 
Uncertainties due to unknown higher-order effects are usually accounted for by vary- 
ing the renormalization (p^) and factorization (pp) scales. In principle, independent 
variations of the two scales should be considered. In fig. ^ we show the scale depen- 
dence of the top cross section. Note that there is compensation of the factorization 
scale dependence when the different subprocesses are added up. In particular, the 
gg and qg factorization scale dependence have opposite behaviour. In fact, it is only 
the combined scale dependence that can be considered an estimate of the neglected 
sub leading corrections. As a second point, we observe that the maximum of the cross 
section is reached around mt/2. Thus, the usual choice of the range mt/2 < p < 2mj 
appears to be particularly justified in this case. As a third point, we observe that the 
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Figure 3: Contribution of Coulomb effects at order and higher, relative 
to the NLO result, for the individual channels and for the total, as a function 
of the CM energy in pp collisions. 



scale dependence in the cross section goes in the same direction for the two scales, 
so that, for the purpose of estimating the associated uncertainty, it is sufficient to 
consider the simultaneous variation of /iR and /ip- 

Aside from the scale uncertainties, which reflect the limitations of the perturba- 
tive QCD calculation, there are uncertainties associated to our imprecise knowledge 
of the physical parameters involved. In particular, the strong coupling constant is 
determined within a certain accuracy. In the case of top production at the Tevatron, 
larger strong couplings tend to give larger partonic cross section. However, for larger 
strong couplings, the Altarelli-Parisi evolution, which occurs along the wide span in 

from the values at which deep-inelastic fits are performed (~ 10 GeV^) to the top 
mass, softens the quark parton densities, thereby decreasing the cross section. In or- 
der to perform a fair estimate of the uncertainty due to Aqcd, we need sets of parton 
densities fitted with different values of the strong coupling. The sets of ref. meet 
our purpose^. In fig. ^ we show the cross section as a function of the strong coupling. 

^The values of A4 that accompany the fortran program for the structure function sets of ref. | |20| 
are not consistent with the values of quoted there, the differences being of the order of 1%. In 
the present work, we extract the values of A5 from their quoted values of as using the standard 
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Figure 4: Scale dependence of the top cross section at the Tevatron. Dotted 
lines are the ji-p dependence at fixed /xr, dashed lines are the /xr dependence 
at fixed fxp, and solid lines are obtained by the simultaneous variation of /xr 
and ■ 
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Figure 5: Top cross section as a function ofas{Mz). The dotted line does 
not include the variation of the parton densities due to the change in Os- 



For comparison we also show the dependence if the parton densities are kept fixed. 
We see the remarkable reduction in the as dependence due to the compensation of 
the rise of the partonic cross section and the decrease of the quark parton densities. 

Our results for the top cross section are collected in tables |I], ^ and for values 
of the top mass that span the current uncertainties. We also show results obtained 
with with the recent parametrization derived in ref. p2| by including the CDF jet 
data. We stress that the numbers in the table are obtained with a standard two-loop 
calculation. No resummation effects are included. 

Our ranges for the top cross section are thus given by 7.78^^'^^ pb, 4.75][^gg2 pb, 
2.97]'^g'3g pb for mt =160, 175, 190 GeV respectively, and are adequately fitted in this 
range by the expression 

— . 175 — JTif I TO 

a{ti)=e^{A.75tl-ll). (3.7) 



As a central value for our determination we have chosen the MRSA' ||2^ result with 
/Ur = /^F = ""^t; in association with a value of A5 = 0.152 GeV (which corresponds to 
as{Mz) = 0.1113, according to the standard two-loop formula [^). For the MRS 
sets with variable A, we have used A5 = 0.0994, .140, .190, .253, .328 GeV (which 

two- loop formula [pl[. 
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CTEQIM 


CTEQ' 


MRSA' 


MRS, variable A, as{Mz) = 


0.105 


0.110 


0.115 


0.120 


0.125 




5.21 


5.07 


5.00 


4.78 


4.99 


5.18 


5.34 


5.48 


mt 


4.96 


4.86 


4.75 


4.57 


4.76 


4.92 


5.05 


5.16 


2mt 


4.38 


4.38 


4.25 


4.13 


4.27 


4.38 


4.47 


4.52 



Table 1: Total cross sections for nit =1'75 GeV at NLO. 





CTEQIM 


CTEQ' 


MRSA' 


MRS, variable A, a,(Mz) = 


0.105 


0.110 


0.115 


0.120 


0.125 


mt /2 


8.64 


8.18 


8.19 


7.80 


8.17 


8.50 


8.81 


9.06 


mt 


8.16 


7.85 


7.78 


7.45 


7.78 


8.06 


8.32 


8.53 


2mt 


7.21 


7.08 


6.95 


6.72 


6.97 


7.18 


7.35 


7.48 



Table 2: Total cross sections for mt =160 GeV at NLO. 
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/^R=/iF 


CTEQIM 


CTEQ' 


MRSA' 


MRS, variable A, a,(Mz) = 


0.105 


0.110 


0.115 


0.120 


0.125 


mt /2 


3.26 


3.22 


3.13 


3.00 


3.13 


3.24 


3.32 


3.39 


mt 


3.08 


3.08 


2.97 


2.87 


2.98 


3.07 


3.14 


3.20 


2mt 


2.72 


2.77 


2.66 


2.59 


2.67 


2.73 


2.78 


2.78 



Table 3: Total cross sections for mt =190 GeV at NLO. 



corresponds to a,(Mz) = 0.105, 0.110, 0.115, 0.120, 0.125). 

In pp collisions at = 2 TeV we obtain, with the same method, a top cross 

pb, 6.53+;;°^ pb, A.utlH 



section of 10.5+!-^ pb, G.^Stif. pb, 4.17+°'^^ pb for mt =160, 175, 190 GeV respec- 



tively. The cross section bands are also shown in fig. |^, for = 1.8 TeV and fig. |^ 
for = 2 TeV. Because of the high precision of this theoretical prediction, the 
measurement of the top cross section could become a sensitive probe of new physics 

PI 



For reference, we also quote the cross section for top production at the LHC We 
get a{ti) = 0.77|'|q^2 ^b, for = 175 GeV and -s/S* = 14 TeV. The error is obtained 
with the same scale and as variations we used for the Tevatron cross section. 
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Figure 6: Top cross section at the Tevatron at ^/S = 1.8 TeV. The solid 
line is obtained with MRSA' parton density, and the dashed lines correspond 
to the upper and lower values obtained in tables 1-3. The experimental data 
are taken from ref. M^. 
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